Bistability and macroscopic quantum coherence in a BEC of 1 Li. 
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We consider a Bose-Einstein condensate of 7 Li in a situation where the density undergoes a 
symmetry breaking in real space. This occurs for a suitable number of condensed atoms in a double 
well-potential, obtained by adding a standing-wave light field to the trap potential. Evidence of 
bistability results from the solution of the Gross-Pitaevskii equation. By second quantization, we 
show that the classical bistable situation leads, in fact, to a macroscopic quantum superposition or 
Schrodinger cat (SC) and evaluate the tunneling rate between the two SC states. The oscillation 
between the two states is called macroscopic quantum coherence (MQC); we study the effects of 
losses on MQC. 



I. INTRODUCTION 

The availability of Bose-Einstein condensates of 
trapped cold atoms ||, has opened the possibility 
of a laboratory engineering of quantum states with a large 
number of atoms Q (around a thousand for 7 Li 0,0]). 

One of the most challenging endeavours of quantum en- 
gineering is the evidence of superposition states [so called 
SC=Schr6dinger cat (SC)] whose mutual interference is 
called macroscopic quantum coherence (MQC). SC have 
been observed, e.g., for states of a trapped ion j?j and of 
a microwave field in a high-Q cavity |g| . 

In this paper we demonstrate the reliable preparation 
of SC consisting of a Bose-Einstein condensate (BEC) of 
7 Li atoms that have negative scattering length, trapped 
in a double-well potential. Realistic calculations have 
been offered for macroscopic quantum tunneling (MQT) 
H [to]. Indeed, combining the kinetic and potential 
terms of a harmonic trap with the inter-particle attrac- 
tion yields a metastable state for N < N c (N c = critical 
population for an attractive BEC). Quantum tunneling 
from this metastable state towards the collapsed state, 
which would otherwise be reached for N > N c , has been 
shown to be feasible. 

A BEC of atoms with negative scattering length, 
trapped in a double well, undergoes a space symme- 
try breaking beyond a threshold number of atoms Ni, 
whereby two stable states are formed. This phenomenon 
has been dealt with theoretically by a two mode ap- 
proach 

We have studied the problem by finding numerically 
the stationary solutions of the Gross-Pitaevskii (GP) 
equation discretized over a space lattice, with reference 
to the 7 Li case |l2|]. Once the stationary solutions have 
been found, we introduce a quantum two mode model, 
with the two modes chosen in such a way as to reproduce 
the stationary solutions of GP. The model shows the fea- 
sibility of macroscopic quantum coherence (MQC). 

Bistability occurs only for an attractive interatomic 



potential (negative scattering length) . Two proposals for 
reaching SC in a BEC with repulsive atoms have been 
put foward |L3[ They both require a Raman cou- 

pling between two different components; Both the pa- 
pers consider copropagating light beams. In Ref. ] ^3[ , 
MQC is shown to require values of the scattering length 
between atoms of different magnetic number substan- 
tially larger than the scattering length between atoms of 
equal magnetic number. This requirement is too strong, 
since no experimental technique is today accessible to 
provide such a difference; furthermore, if such a differ- 
ence could be achieved, outstanding symmetry-breaking 
effects would occur . Reference [Q introduces a time- 
dependent evolution, so that SC is reached over a time 
of the order of 1 sec. However, no clear-cut experimental 
test is offered to discriminate between SC and a statisti- 
cal mixture of two separate states. 



II. BISTABILITY AND SYMMETRY 
BREAKING 

We refer to a condensate of 7 Li atoms trapped in a 
double-well potential. A suitable model for it is given by 



V{x) = -m 



2 2, 2/2 
LOnX 1 + U!j_(X 2 



■ A-cos(2tt— ) (1) 



The quadratic part is due to the interaction of the atoms 
with the magnetic field of the trap. According to labora- 
tory implementations we choose u>\\ = 2tt x 130s -1 and 
u>± = 2ir x 150s -1 . The additional term is generated by 
two opposite laser beams in a standing wave configura- 
tion. A suitable choice of the standing-wave parameters 
yields a double well-potential. Taking into account the 
interatomic interaction, the atomic system is described 
by a macroscopic wave function ijj that satisfies the GP 
equation 



ih 



~dt 



2m 



y^ + ViJ; + g\^=(H + g\^\ 2 )^, (2) 



2 



with H = —h 2 /(2m)\/ 2 + V. Here, m ~ 7a. u. is the mass 
of the lithium atom, and g — {Ai:Ti 2 /m)a s , where a s is 
the s-wave scattering length for 7 Li, a s ~ — lAbnm. For 
a small number of atoms the GP nonlinearity can be ne- 
glected and Eq. (||) reduces to an ordinary Schrodinger 
equation. In such a case and for a sufficiently high bar- 
rier, the lowest energy level is described by a two-peak 
wave function symmetric with respect to inversion of the 
space axes, once the coordinate origin coincides with the 
trap center (Fig. [I| dashed-dot line). 

Figure [l] reports the spatial distribution of the ground 
state of Eq. (||) for different numbers of trapped atoms; 
we have used a numerical method that consists in solving 
GP on a discrete space lattice and evaluating the lowest 
energy state. The barrier is specified by the two param- 
eters A n = A/h — 2650s -1 and a = 5/j,m (see Eq. (1)). 
As shown in the figure, for N = 450 the distribution is 
symmetric; instead for N = 500 the nonlinear term is 
sufficient to destabilize the symmetric state, giving rise 
to two asymmetric stable states. For N — 600 one well is 
almost empty. As we increase the number of atoms, the 
nonlinearity plays a relevant role. By a self-consistent 
argument we realize that the symmetric wave function 
becomes unstable and we can have two new minimal en- 
ergy states with distribution no longer symmetric for in- 
version (symmetry breaking). Indeed, let us assume a 
distribution as in Fig. [j] (dashed or solid line); then the 
effective potential for such a distribution, due to the sum 
of the external potential with g\ip\ 2 is an asymmetric dou- 
ble well with the lower minimum corresponding to the 
higher population peak. For a sufficiently high nonlinear 
term the potential imbalance stabilizes the asymmetric 
distribution as in Fig. [jj 

We confirm the numerical calculation by the following 
analytic model. Let ip a be the equilibrium symmetri- 
cal wave function (either stable or unstable), and ip b be 
a suitable antisymmetrical wave function such that the 
weighted sum of the two wave functions lowers either one 
of the two peaks. In the two-dimensional space of these 
wave functions, any other one can be expressed as 



ip(Hi) = aip a (x) + btpb(x). 



(3) 



Without loss of generality we can choose ip a and 4>b 
as real functions and thus consider real values for a and 
b. As we take / il> 2 d 3 x — J tp b 2 d 3 x = 1, it follows that 
a 2 + b 2 = Jtp 2 d 3 x = M(a,b), where M(a,b), the total 
number of atoms in the condensate, depends upon a and 
b. The energy is then given by 

H = a 2 H aa + b 2 H bb + \{a A I aa + b A I bb ) + 3a 2 b 2 I ab (4) 



where 



r a H^ a d 3 x, H, 



bb 



I 



^Htfj b d 3 x, 



laa = 9 / ^ a d xx, I bb = g tj} h d X, 



iy(x,o,o)i 
do-'V 3 ) 

1.5 



0.5 




N=450 
N=500 
N=600 



FIG. 1: Density distribution of the lithium BEC in real space 
for N = 450 (dash-dotted line), N = 500 (dashed line), and 
N = 600 (solid line) with A n = A/h = 2650s -1 and a = 5(j,m. 
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We look for the minimal of energy with the constraint of 
a fixed number N of condensed atoms. These constrained 
minima satisfy the relations 



d(H - fiM) _ d(H - fiM) 



da 



db 



= 0, 



(5) 



where /j, is a Lagrange multiplier. We thus solve for a and 
6, with the condition N — a 2 + b 2 . The solutions a = 
and 6 = correspond, respectively, to the antisymmetric 
function and the symmetric one. The other solutions, for 
both a and b nonzero, yield a 2 , b 2 values as functions of 
fi. Using the constraint of fixed TV, we eliminate \i and 
find 



H, 



bb 



N{3I ab 



ibb 



b 2 = 



Slab — laa — Ibb 
H bb ~ H aa + N(3Ig b - lag) 
Slab — laa — Ibb 



(6) 
(7) 



Since g is negative, the denominator is always negative. 
Indeed, ip 2 and ip 2 are almost equal at each point x. For 
low N, the dominant terms in the numerators have oppo- 
site sign, thus one of the two squares has to be negative, 
which means that there is no solution with asymmetric 
wave function. On the other hand, for sufficiently large 
N, the second term in the numerator of the two equations 
can compensate for the positive one, since the quantities 
3I ab — I aa , 3I ab — I b b are always negative. Thus, we have 
proved that two asymmetrical steady states exist beyond 
a threshold value of N. 

In order to prove that the two states are stable, it is 
sufficient to show that the symmetrical state becomes 



3 



unstable above threshold, that is, 



d 2 H 



db 2 



< 



(8) 



b=0 



where we have a 2 = N — b 2 . It is easily found that 



d 2 H 



db 2 



2{H bb - H aa ) + 2N('SI ab - I aa ) (9) 



b=0 



Going back to Eq. (^) , the two stationary asymmetrical 
states occur when the numerator changes sign. In fact, 
the numerator is the right side of Eq. (||), hence, at the 
critical point the symmetrical state becomes unstable. 

One can easily evaluate the threshold value Ni of N 
for which this symmetry breaking occurs. 



III. MACROSCOPIC QUANTUM COHERENCE 

In such a bistable situation the energy displays two 
minima of equal value in the infinite-dimensional phase 
space. From a classical point of view, as the system is 
in its lowest energy state, the condensate is localized in 
either one of the two minima, where it will remain in the 
absence of thermal noise once we keep the atom number 
constant. Since, however, the condensate is a mesoscopic 
system, quantum fluctuations play a relevant role. This 
can be shown by second quantization of the field, re- 
placing the c-number macroscopic wave function by field 
operators. 

Quantum fluctuations allow the passage from one to 
the other stable state without thermal activation, by pure 
quantum tunneling. Furthermore, due to the coherent 
nature of the process, we expect coherent oscillations be- 
tween the two wells, that is, MQC. We now evaluate the 
tunneling rate as a function of the system parameters 
showing the feasibility of MQC. 

The most natural way of evaluating the tunneling 
rate consists in finding the two lowest eigenvalues of the 
Hamiltonian and taking their difference. Indeed, the sum 
and difference of the corresponding states are respectively 
the alive and dead states of SC, and the transition time is 
half the period corresponding to the frequency difference. 

The problem is simplified by reducing it to two degrees 
of freedom by the expansion of Eq.(||). 

We select the basis functions tp a and ip b as follows. 
Calling ipo(x) the minimal-energy wave function of the 
CP problem (see e.g. Fig. [l]), we take for ip a and 
ipb, respectively, the symmetric and antisymmetric sums 
ipo(x) ±tpo(—x). Expansion (|J) with these ip a and ip b in- 
cludes the original functions ipo(x) and ipo(~ %) f° r suit- 
able values of a and b. Furthermore, it simplifies the form 
of the Hamiltonian, as we see right now flq| . 

In second quantization, a and b in Eq. (3) become the 
operators a and 6, obeying Bose commutation rules with 



(a) 



- N=621; A n =2650 s" , c=4.5 |im 
■ N=496; A =2650 s"', c=5 flm 




- N=765; A n =2000 s" , c=4.5 flm 

- N=642; A =2000 s"', o=5 flm 



FIG. 2: First two excited energy levels versus the number 
N — Ni of atoms above the threshold Ni. Within each figure 
we keep the laser amplitude fixed and just vary the pitch of 
the potential lattice; (a) corresponds to A n = A/Ti — 2650s -1 
and (b) corresponds to A n = 2000s" 1 . 



their conjugates d\ w . Exploiting the operator version 
of Eq. (3) and its adjoint, the Hamiltonian becomes 

H = a)aH aa + tfbH bb + i(a t2 a 2 / QQ + P 2 b 2 I bb )+ (10) 



\^ 2 b 2 + i& t2 a 2 + 2a)tfdb ) I ab 



where the coefficients H aa , H bb , I aa , I bb and I ab are 
the same as in Eq. (4). 

We consider the basis of eigenvectors of the number 
operators 



|0,iV>,|l,JV-l >,...|JV,0> 



(11) 



where a^a|fc, m >= k\k, m > and b*b\k, m >= m\k, m >. 
Let us call 



Hi M =< l,N -l\H\k,N -k > 



(12) 



the generic matrix element of the Hamiltonian on the 
above basis. We evaluate the eigenvalues of this matrix. 
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FIG. 3: Same as with Fig. ^[ but with fixed a, respectively 
(a) a = 4.5/im and (b) a = 5/im. 



We have considered two different wavelengths and two 
amplitude values of the applied field. 

In Fig. H we report on the first two excited energy lev- 
els versus the number of condensed atoms beyond the 
threshold value. In such figure we keep constant the bar- 
rier height and just vary the barrier width. Notice that 
for increasing a the maximum tunneling frequency re- 
duces, but the slope at which it reduces for increasing 
N — Ni is less steep. 

In Fig. ^ we keep fixed the barrier width and change 
its height. Here too the maximum frequency decreases 
for increased heights, but again as in Fig. g, the slope 
decreases for increasing N — Ni. 

In Fig. [| we reduce the threshold for BEC breaking 
by increasing lu±. Precisely, we report on the first two 
excited states for ui± = 2tt x 600s -1 , A n = 2000s~ 1 and 
a = 5/im. The threshold value is N ~ 190. 

In the number representation there is no explicit ev- 
idence of a SC as a two-peak distribution (see Fig. ||). 
We look for a suitable observable, whose probability 
distribution provides such an evidence. We take the 
first component of the barycenter coordinate x c — 
(1/N) J xi\ip\ 2 d 3 x as the appropriate variable since the 
corresponding classical states [minima of the Hamilto- 



FIG. 4: First two excited energy levels versus the number of 
atoms N. The system parameters are wj_ = 2ty x 600s -1 , 
W|| = 2tt x 130s"" 1 , An = A/h = 2000s -1 and a = 5/im. 
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FIG. 5: Probability distribution of the population in the 
ground state as a function of the number of atoms in mode 
'b'. A = 2000s -1 , cr = 5/im and N = 655. 



nian (^)] have separated barycenters. It is associated 
with the operator 



1 

N 



x 1 ^ a (x)ip b (x)d 3 x(a^b + atf) (13) 



Besides a c-number factor, the observable is thus M = 
a)b + aW. The associated probability density is P(m) = 
< to | </>o > | 2 , where \(j> > is the ground state of Hamil- 
tonian ( |io| ) and |to > is the eigenstate of M with eigen- 
value to. Notice that once the number of atoms is fixed 
the eigenvectors are not degenerate. 

Since we know the components of \(f>o > on the number 
basis we must express |m > with respect to that basis, 
that is, 



to >= 



E 

k 



cT\k,N~k> 



(14) 
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FIG. 6: Probability distribution of the "barycenter" indicator 
for the lithium BEC, for different N values and for a standing 
wave with A n = 2000s -1 and a — 5fj,rn; (a) lj± = 2ttx 150s -1 
and (b) u± = 2ir x 600s -1 . 



Applying the operator M to the above ket we have 

m\m>=^ c T(a r b + ab r )\k,N - k> . (15) 
k 

Projecting on < I, N — l\ the above ket it follows that 

Y,Mi, k c% = mc?, (16) 



where 



M hk =< l,N — l\M\k, N-k> 



(17) 



Since Mi jk are known, solving the eigenvector equation 
( |l6| ) we can evaluate the coefficients c™ and hence P(m) . 

The two-peak distribution P(m) is plotted in Fig. g(a) 
for different N values and for A/K= 2000s -1 , a = 5/xto. 



The same distributions, but for u>± = 2tt x 600s -1 , are 
given in Fig. ^(b). 

Experimentally, evidence of SC against a trivial statis- 
tical mixture is obtainable by setting the system at the 
above parameter values and observing the coherent oscil- 
lation between the two energy minima by a measurement 
not destroying the coherence (e.g., a phase contrast tech- 
nique). The condensate can be prepared in the ground 
state, because the system condensates naturally in this 
state. Then one measures the condensate barycenter by 
the above noninvasive technique at the initial time. This 
measurement collapses the distribution P(m) onto one 
peak. If the energy transfer of the measurement is not too 
large, we expect that only the first excited state is popu- 
lated and hence we have the superposition of <f>a(m) and 
0i(to) [see Fig. 0(a)] . At the following times the distri- 
bution oscillates between the two peaks. The symmetry 
of Hamiltonian ( |l0| ) with respect to the transformation 
a — * —a shows that 4>o{rn) and 0i(rn) are symmetric or 
antisymmetric; by a suitable choice of the phase, both 
0o and 0i can be real functions. We consider a state 
preparation such that at t — the state 4>q (to) + <f>i (m) 
has only one probability peak [see Fig. 0(a), dashed line]. 
After a quarter period corresponding to the frequency 
separation cj/ between ground and first excited states, 
the superposition will be 0o (to) — i(f>i (to) and the corre- 
sponding probability is the two-peak aspect (solid line), 



Q(m,t = n/ujf) = 4>o{m) 2 + 4>i(mY 



(18) 



At time t = n/u>f the only peak is that missing at time 
t = 0, thus there is a coherent oscillation between the 
two states. Detecting such an oscillation would provide 
evidence of an SC at an intermediate time when both 
peaks are present. 

The intermediate distribution, given by Eq. (|l8|), may 
display two peaks even when the ground state has only 
one [see Fig. |(b)]. This occurs because the distribution 
0i (to) 2 has always two peaks. 

Notice that if we had N loosely coupled or independent 
atoms (Ng sufficiently low or even zero) the superposition 
of ground and first excited state would have a single peak, 
oscillating as a coherent state inside a harmonic poten- 
tial. This would by no means be a SC. On the contrary, 
we have shown that for Ng sufficiently high, we have a 
two peak distribution with the two partial barycenters at 
nearly constant positions. During the evolution, the two 
peak amplitudes oscillate, that is, the probability to find 
the system in either state oscillates. 

In conclusion, we have found that the oscillation fre- 
quency between the two states of the SC is u/f ~ 50s -1 
sec Fig. |(a) for N—Ni - 5]. In order to neglect thermal 
activation, the system should be cooled at a temperature 
of around InK. To cool at InK is within the reach of 
present technologies, even though such a low tempera- 
ture has not yet been reported. On the other hand, the 
collective degrees of freedom for which there is quantum 
coherence may be weakly coupled with the other modes 
of the condensate, which act as a thermal bath. Hence, 
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FIG. 7: Time evolution of the probability distribution. Wf is 
the tunneling frequency. The parameters of the standy wave 
are as in Fig. [](a); (a) number of atoms N = 655 and (b) 
N — 645, for which the ground state does not display two 
separate peaks [see solid line of Fig. 0(a)] . 



in such a case it might not be necessary to cool at InK 
the whole condensate, but just the involved degrees of 
freedom. 



IV. LOSSES AND OBSERVABILITY OF MQC 

Besides the finite temperatures effects, we must also 
account for atom losses from the condensate. 

The entanglement between condensate and lost atoms 
will imply a decoherence of the superposition. Near the 
threshold for symmetry breaking, the two macroscopic 
wave functions corresponding to the energy minima of 
the classical system are almost coincident (see Fig. [l], for 
N = 500), thus the loss of one or a few atoms does not 
allow us yet to discriminate between the two alternatives, 
and thus decoherence effects are almost negligible. 

We perform a simple calculation of the decoherence 
effect due to atom loss for the condensate. In the Hartree- 
Fock approximation all the atoms are in the same state, 
given by the GP wave function. Hence, if tpo(x) is one 



of the two asymmetrical states with minimal energy, the 
wavefunction for the whole system of N atoms is given 
by 



JY 



ipi(xi,x 2 , ...,x N ) = JJVo(^i)- 



(19) 



But the GP yields another minimum ipa(~ x), for which 
the overall wave function is given by 



1pr(xi,X 2 , ...,X N ) = Y]_ i-'o(-Xi). 



(20) 



The coherent superposition of the two alternatives is 
given by 



lp s = K(lpl + 1p r ) 



(21) 



where K is a normalization factor. If ?p^ r are orthogonal 
then K = -4=. 

v2 



The density operator is 



1 



p= o(l^) + hM)((^| + (Vv|) 



(22) 



The coherence between the two alternatives is given by 

C = 2Tr[|^>(Vv|p] = l. (23) 

If we assume that the atoms escaping the condensate 
transfer a negligible energy to the trapped atoms, then 
the residual coherence after the loss of M atoms is 



C = 2Tr[$ l )(i> r \p} ) 



(24) 



where the vectors \ipi, r > refer to the N — M atoms still 
in the condensate. The functions ij)i r can be written as 



M 



Lr 



4>l,r JJ^o(±^i), 



(25) 



where i = 1, 2, M are the lost atoms. We rewrite p by 
using this expansion, and find, provided that (ipi\tpr) = 0, 



C = 



(Fxip (x)il> (-x) = (I ) 



A I 



(26) 



Let us call e = 1 — I . If e <C 1 then the coherence is 
given by 



C 



-eM 



(27) 



The quantity N e i = 1/e provides the number of atoms 
that must be lost in order to reduce the coherence by 1/e. 
Figure || shows how Nd scales with the total number of 
atoms N for A = 2000s -1 and a = 5pm. 

Notice that near the threshold value (N — 643) Nd is 
rather large, and hence, coherence is more robust with 
respect to atom losses. Far threshold Nd has a weak 
dependence on N. 
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FIG. 8: Number Nd of lost atoms for which the coherence 
decays by 1/e versus the total number of atoms, for A„ = 
A/h = 2000s" 1 and a = 5fim. 



FIG. 10: Number Nd of lost atoms for which the coherence 
decays by 1/e versus the total number of atoms, for the pa- 
rameters of Fig. 0. 
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FIG. 9: R(N) for A n = A/h = 2000s" 1 and a = 5/j.m. 



We now evaluate the loss rate of the condensate. The 
relevant processes are two-body inelastic and three-body 
collisional decays. The total loss rate of the two decays 



are [17 



R(N) = aN^ J c£V|^ (r)| 4 + LN A J d A r\^{x)\ b (28) 

where a is the two-body dipolar loss rate coefficient and 
L is the three-body recombination loss rate coefficient. 
We use the following values: a — 1.2 x 10" 14 cm 3 s -1 [fj"7| , 
L = 2.6 x 10 -28 cm 6 s -1 pj. 

We report in Fig. | R{N) for A = 2000s" 1 and a = 
5fim. From Figs. ^(b)|| and [| we can infer that near 
threshold atom losses have a small effect on coherence. 

For instance, for N = 650 the tunneling rate is about 



25s" 1 [see Figs. ||(b) and||(b)], thus the number of atoms 
lost during an oscillation period (~ 0.25s" 1 ) is around 
3.5 (see Fig. ^); but Fig. || shows that such a loss is 
negligible for coherence. For N = 655 the tunneling rate 
is about 10s _1 , thus the number of atoms lost during an 
oscillation period is about 9, still below ~ 30. 

Besides decoherence, the atom loss yields a shift in 
the tunneling frequencies. As shown in Figs. ^-||, the fre- 
quencies vary a lot as we reduce the number of condensed 
atoms. For N — 650 the lost atoms are around 3.5 dur- 
ing an oscillation period, thus the oscillation frequency 
changes by ~ 30%. 

In order to reduce the effect of losses we lower the 
threshold value by increasing the frequency ujj_. In 
Figs, [h] and [ll] we report on Nd and R{N) for u± = 
2ir x 600s -1 ; the longitudinal frequency and the stand- 
ing wave have been kept unchanged. The tunneling fre- 
quencies and probabilities P(m) are given in Figs. ||(b) 
and |6](b) , respectively. We notice that here the loss rate 
is much smaller, and hence the frequency shift is smaller, 
whereas the tunneling frequencies are slightly larger. A 
further reduction of the threshold value should improve 
the situation. 



V. CONCLUSION 

We have solved numerically the Gross-Pitaevskii equa- 
tion for a 7 Li condensate in a double-well potential and 
we have shown that the spatial density undergoes a 
symmetry breaking for a suitable number of condensed 
atoms. The classical asymmetrical wave function is used 
to select two modes. With their quantization we have 
evaluated the tunneling rate for some parameters of the 
system and we have found that it is within the reach of 
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FIG. 11: R(N) for the parameters of Fig. §. 



present laboratory technologies. Thermal effects could be 
lower than expected because of a decoupling between the 
collective variables and the thermal bath. To prove the 
presence of two alternatives we have introduced an ap- 
propriate observable. We have then evaluated the effects 
of loss and we have showed that the decoherence is neg- 
ligible. In a forthcoming work we discuss the conditions 
for bistability in the case of mutual repulsion (positive 
scattering length). 



[1] 
[2] 

[3] 
[4] 

[5] 
[6] 



[7] 
[8] 
[9] 

[10] 

[11] 
[12] 



M. H. Anderson et al., Science 269, 198 (1995). 
K. B. Davis et al., Phys. Rev. Lett. 75, 3969 (1995). 
C. C. Bradley et al., Phys. Rev. Lett. 75, 1687 (1995). 
C. C. Bradley, C. A. Sackett, and R. G. Hulet, Phys. 
Rev. Lett. 78, 985 (1997). 

R. Walsworth and L. You, Phys. Rev. A 56, 555 (1997). 

In the case of positive scattering length as for 87 Rb Jl|] or 

23 N a j^], the number of atoms in a condensate may be 

as large as 10 6 or 10 7 , respectively. 

C. Monroe et al., Science 272, 1131 (1996). 

M. Brune et al., Phys. Rev. Lett. 77, 4887 (1996). 

M. Ueda and A. J. Leggett, Phys. Rev. Lett. 80, 1576 

(1998). 

L. Salasnich, Phys. Rev. A 61, 15601 (1999). 



Tin-Lun Ho, and C. V. Ciobanu, 
In a previous report, "A. Montina and F. T. Arecchi, J. 



:ond-mat/0011095 



Mod. Opt., 49, 319 (2002)", we had considered this prob- 
lem but in less detail; furthermore, we did not consider 
the role of losses on the observability of MQC. 

[13] J. I. Cirac et al., Phys. Rev. A 57, 1208 (1998). 

[14] D. Gordon and C. M. Savage, Phys. Rev. A 59, 4623 
(1999). 

[15] D. Gordon and C. M. Savage, Phys. Rev. A 58, 1440 
(1998). 

[16] A generic two mode model is also used by C. Menotti 
et al., Phys. Rev. A 63, 023601 (2001), however, here 
we have selected the suitable modes by reference to the 
associated GP equation. 

[17] R. J. Dodd et al., Phys. Rev. A 54, 661 (1996). 

[18] A. J. Moerdijk et al., Phys. Rev. A 53, 916 (1996). 



